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Abstract 

We give a new proof of the KAM theorem for analytic Hamiltonians. The proof is 
inspired by a quantum field theory formulation of the problem and is based on a 
renormalization group argument treating the small denominators inductively scale 
by scale. The crucial cancellations of resonances are shown to follow from the Ward 
identities expressing the translation invariance of the corresponding field theory. 



1 Introduction 



Consider the Hamiltonian 

H(I, (f>) = u-I+^I-fil + X U(cf>, I) (1.1) 

with <f> G R d /(27rZ d ) = T d , I G R d , u G H d with the components Ui independent over 
Z d and \x a real symmetric d x d matrix. It generates the Hamiltonian flow given by the 
equations of motion 

= cj + /i/ + AcW, I = -Xd^U. (1.2) 
l Partially supported by EC grant CHRX-CT93-0411 
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For the parameter A = and the initial condition (0o, 0), the flow (4>Q + ut, 0) is quasiperi- 
odic and spans a <i-dimensional torus in T d x R d . KAM-theorem deals with the question 
under what conditions such quasiperiodic solutions persist as the parameter A is turned 
on. 

Consider a quasiperiodic solution in the form 

((j)(t),I(t)) = (00 + ^ + 9(00 + ^), J(0o + W*)). 

Eqs. ( |OD require that Z = (6, J) : T d — > R d x R d satisfies the relation 

VZ(<f>) = -A dU(<f> + 0(0), J(0)) , (1.3) 



2>=( °* (1-4) 



where 9 = (<9^, 9/) and 

-a; • ^ /i 

Note that if Z is a solution of Eq. (|1.3| ) then so is Zp for /3 e R d and 

^(0) = Z(0-/3)-(/3,O). (1.5) 

Eq. (|1.3|) is a fixed point problem for the function Z of a difficult type: the straightforward 
linearization D + XddU(<j),0) is not invertible for any interesting [/ (see e.g. ||). Also, 
one can expect to have a solution only for sufficiently irrational uo e R d , e.g. satisfying a 
Diophantine condition 

|w • q\ > a\q\~ v for qeZ d ,q^0 (1.6) 

with some g-independent a, v > 0. There have been traditionally two approaches to the 
problem: 

1. The KAM approach. (|1.3|) is solved by a Newton method that constructs a sequence 
of symplectic changes of coordinates defined on shrinking domains that, in the limit, 
transform the problem to the A = case [[j], 0, |17|, [18| . 

2. Perturbation theory. For U analytic (see below) one can attempt to solve ( |1.3| ) by 
iteration. This leads to a power series in A, the Lindstedt series: Z = J2 n Z n X n . Each Z n is 
given as a sum of several terms (see Sect. 9), some of which are very large, proportional to 
(n!) a with a > 0, due to piling up of "small denominators" (w-q) from the momentum space 
representation of operator T)~ x . However, the KAM method also yields the analyticity 
of Z in A ||19[| . Thus the Lindstedt series must converge. To see this directly turned 
out to be rather hard and was finally done by Eliasson [Q who, by regrouping terms, 
was able to produce an absolutely convergent series that gives the quasiperiodic solution. 
Subsequently Eliasson's work was simplified and extended by Gallavotti 
Chierchia and Falcolini R U\ and by Bonetto, Gentile, Mastropietro |12|, |13, 




1ft by 
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In the present paper we shall develop a new iterative scheme to solve Eq. ( |1.3| ). It is 
based on a direct application of the renormalization group (RG) idea of quantum field 
theory (QFT) to the problem. The idea is to split the operator T> (or rather its inverse, 
see Sect. 2) into a small denominator and large denominator part, where small and large 
are defined with respect to a scale of order unity. The next step is to solve the large 
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denominator problem which results in a new effective equation of the type ( |1.3j ) for the 
small denominator part, with a new right hand side. The procedure is iterated, with the 
scale separating small and large at the n th step equal to rf 1 for some fixed rj < 1. As a 
result we get a sequence of effective problems that converge to a trivial one as n — > oo. A 
generic step is solved by a simple application of the Banach Fixed Point Theorem in a big 
space of functionals of Z representing the right hand side of Eq. ( |1.3|) in the n th iteration 
step. 

Our iteration can be viewed as an iterative resummation of the Lindstedt series, as 
will be discussed in Sect. 9. This iterative approach trivializes the rather formidable 
combinatorics of the small denominators. The functional formulation in terms of effective 
problems removes also the mystery behind the subtle cancellations in the Lindstedt series: 
they turn out to be an easy consequence of a symmetry in the problem as formulated in 
terms of the so called Ward identities of QFT. The QFT analogy of the problem (|1.3|) has 
been forcefully emphasized by Gallavotti et al. [jllj, [12]]. The proof of Eliasson's theorem 
by these authors was based on a separation into scales of the graphical expressions entering 
the Lindstedt series and was a direct inspiration for the present work. 

An important part of the standard RG theory is an approximate scale invariance of 
the problem that is exhibited and exploited by the RG method. The KAM problem also 
is expected to have this aspect: as the coupling A is increased the solution with a given 
oj eventually ceases to exist. For suitable "scale invariant" oj (e.g. in d = 2 for oj = (1,7) 
with 7 a "noble" irrational) the solution at the critical A is expected to exhibit a power law 
decay of Fourier coefficients and periodic orbits converging to it have peculiar "universal" 



scaling properties [16 , 21]. We hope that the present approach will shed some light on 



these problems in the future. 

While the main goal of this paper is to develop a new method, we use it to reprove 
the following (classical) result: 

Theorem 1. Let U be real analytic in <p and analytic in I in a neighborhood of I = 0. 
Assume that oj satisfies condition (\1.6j ). Then Eq. ( |1.3|) has a solution which is analytic 
in A and real analytic in <ft provided that either 

(a) (the non-isochronous case) /i is an invertible matrix and |A| is small enough (in a 
^-dependent way). 

(b) (the isochronous case) [l = 0, / d diU(<p,0) dip = 0, the dx d matrix with elements 
J T d d^diJJi^cf), 0) d(j), k, I = 1, . . . , d, is invertible and |A| is small enough. 

The above solutions are unique up to translations j \l.dj ). 

Remark. Actually, we show that the solution is an analytic function not only of A, but 
of the potential U, when the latter belongs to a small ball in a Banach space of analytic 
functions (see Sect. 3 for the introduction of such spaces). This allows us to consider more 
general Hamiltonians of the form 

if (7,0) = 77 (7)+ U(<j>,I). 

with 77 and U analytic and U small. Indeed, we may expand 77 around 7 s.t. <9/77 (7 ) = 
oj, with oj satisfying condition ( |L~6|) : 

77 (7) = 77(7 ) + oj ■ (7 - 7 ) + \ (7 - 7 ) • /x(7 - 7 ) + 77 (7) 
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and define U = U + Hq so as to include in it all the terms of order higher than two in 
the expansion of Hq. Replacing / — 1$ by /, we may apply Theorem 1 provided that U 
satisfies the corresponding hypotheses. Also, more general cases where \i is a degenerate 
matrix can be treated. 

The organization of the paper is as follows. In Sect. 2 we explain the RG formalism. 
In Sect. 3, we introduce spaces of analytic functions on Banach spaces; such spaces will 
be used to solve our RG equations. In Sect. 4, we state the main inductive estimates 
which are proved in Sect. 6 after an interlude on the Ward identities in Sect. 5. Theorem 
1 is proved then in Sect. 7. Sect. 8 explains the connection of our formalism to QFT for 
those familiar with the latter. We should emphasize that the QFT is solely a source of 
intuition, the simple RG formalism of Sect. 2 is independent of it. Finally, in Sect. 9, the 
connection with the Lindstedt series is explained. 



2 Renormalization group scheme 

In this section we explain the iterative RG scheme without spelling out the technical 
assumptions that are needed to carry it out. We refer the reader to Sect. 9 for a graphical 
representation of the main quantities introduced here. 

We shall work with Fourier transforms, denoting by lower case letters the Fourier 
transforms of functions of <f), the latter being denoted by capital letters: 

F{4>) = £ e^ f(q) , where f(q) = [ e 1 ^ F{cj>) dcfi 
with dcf) standing for the normalized Lebesgue measure on T d . 



Note first that we may use the translations ( |1.5| ) to limit our search for the solution 
of Eq. ( |1.3| ) to the subspace of with zero average, i.e. with #(0) = in the Fourier 
language. It will be convenient to separate the constant mode of J explicitly by writing 
Z = X + (0, Q where X has zero average. Let us define 

W (<f>; X, C) = A dU (((/>, + X{<f>)) . (2.1) 

Denote by Go the operator — acting on R M -valued functions on T d with zero average. 
In terms of the Fourier transforms, 



H{u-q) 2 i(wq) 



-i 



for q ^ and (G x)(0) = 0. Writing Eq. (|1.3|) separately for the averages (i.e. q = 0) 
and the rest, we may rewrite it as the fixed point equations 

X = G PW (X,C), (2.3) 

(0,/iC) = - [ WoifrX^OdtP, (2.4) 
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where P projects out the constants: PF = F — J F((p)d<j). Our strategy is to solve 

Eq. ( |2.3| ) by an inductive RG method for given (. This turns out to be possible quite 
generally without any nondegeneracy assumptions on U. The latter enter only in the 
solution of Eq. Q2.4|). Below, we shall treat W given by Eq. ( |2.1| ) as a map on a space of 
R M -valued functions X on T d with arbitrary averages^. The vector ( will be treated as 
a parameter and we shall often suppress it in the notation for Wq. 



For the inductive construction of the solution of Eq. ( |2.3| ), we shall decompose 

G = Gi + T , (2.5) 

where Tq will effectively involve only the Fourier components with \u ■ q\ larger than 0(1) 
and Gi the ones with \u ■ q\ smaller than that (see Sect. 4). In particular, we shall have 
T = r -P. Upon writing X = Y + Y, Eq. ( [2.3|) becomes 



Y + Y = (d + ^PWoiY + Y). (2.6) 

Suppose that Y = Y where Y solves for fixed Y the "large denominator" equation: 

Y = F W (Y + Y ). (2.7) 

Then Eq. ( |2.6|) reduces to the relation 

Y = G x PWx{Y) (2.8) 

if we define Wi(Y) = Wq(Y + Yq). We have thus reduced the orginal problem ( |2.3|) to the 
one from which the largest denominators were eliminated, at the cost of solving the easy 
large denominator problem (2.7) and of replacing the map Wq by W\. 

Note that, with these definitions, Y = roWi(y) and thus W\ satisfies the fixed point 
equation 

Wi(Y) = WoiY + ToW^Y)). (2.9) 

Conversely, this equation, which we shall solve for W% by the Banach Fixed Point Theorem 
in a suitable space, implies that Yq = TqWi(Y) satisfies Eq. (\2.7j ) and thus that 



X = Y + T W 1 (Y) = F t (Y) (2.10) 
is a solution of Eq. ( |2.3| ) if and only if Y solves Eq. (|2.8|). 

After 7i — l inductive steps, the solution of Eq. ( |2.3|) will be given as 

X = F n -i(Y), (2.11) 

where Y solves the equation 

Y = G n _iPW n _i(Y) (2.12) 



'That the solution X of Eq. (2.3) has zero average follows from the form of the equation. 
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and G n -i contains only the denominators \u ■ q\ < 0(r] n ) where rj is a positive number 
smaller than 1 fixed once for all. The next inductive step consists of decomposing G n _i = 
G n + T n _i where T n _i involves \u ■ q\ of order if 1 and G n the ones smaller than that. We 
define W n (Y) as the solution of the fixed point equation 

W n (Y) = Wn.^Y + T n ^W n (Y)) (2.13) 

and set 



F n {Y) = F n _ x (y + r n . 1 w n (Y)) 



(2.14) 



(which is consistent with relation ( |2.10|) if we take Fq(Y) = Y). Then replacing Y in 
Eqs. ([HI]) and (gig ) by Y + T^W^Y), we infer that X = F n (Y) if Y = G n PW n (Y) 
completing the next inductive step. Note also the cumulative formulas that follow easily 
by induction: 



W n (Y) = W (Y + T <n W n {Y)) 
F n (Y) = Y + T <n W n (Y) , 



(2.15) 
(2.16) 



ra-l 



where T <n = Tk contains all the denominators larger than 0(r] n ). 

k=0 



We shall control the transformations (|2.13 ) and ( p,14|) in suitable norms. The point 
of the inductive procedure is that PW n (Y) becomes effectively linear in Y for large n, see 
Remark 1 after Proposition 3 below, so that Y = is a better and better approximation 
to a solution of the equation Y = G n PW n (Y). In fact, as follows from the cumulative 
relations (|2.15| ) and (|2.16|) , X n = F n (0) = T <n W n (0) solves the approximate problem: 



X n — T <n W (X„ 



(2.17) 



obtained from Eq. 



by replacing G by T <n (since T <n P = T <n ). We shall construct 



the solution X of Eq. ( |2.3j ) as the limit of the approximate solutions 



X 



lim X n 

n— >oo 



(2.18) 



3 Spaces 

Let us rewrite the definition ( [2.1|) in terms of the Fourier transforms: 

w (q; y) = X I e^ dU((<j>, C) + F(0)) d<j>, 

where we recall that y refers to the Fourier transform of Y. Let us explain here how to 
view Wo as an analytic functional on a suitable Banach space. The analyticity of U implies 
the following. There exist p > 0, a > and b < oo such that the coefficients U m+ i((f), Q, 
belonging to the space of m-linear maps C(C 2d , . . . , C 2d ; C 2d ), of the Taylor expansion 

oo 

dU{{cj>, C) + Y) = J2 i£W</>, CW, ■ ■ ■ , Y) 

m=0 
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are analytic in |£| < p and their Fourier transforms satisfy the bounds 

E^ kl IK+i(<7,C)ll £(c2d c2d;C2d) < bm\p- m . (3.1) 

q 

For later convenience, we shall use in C 2d = C d x C d the norm | ■ | defined by 

\{zx,z 2 )\ = \zi\ + \z 2 \ (3.2) 

and the induced norms on the spaces of linear maps. Inserting the Fourier series for Y 
we end up with the expansion 

oo 

w (q;y) = E E^ M ™+i(9 - E*>0(y(<?i),---,2/(<?m)) 

m=0 q 

oo 

= E E U, m) (9,9l,---,9m;C)(l/(9l),---,Z/(9m)), (3-3) 
m=0 q 

where q = (qi, ■ ■ ■ ,q m ) G Z md . This formula suggests to consider w as an analytic 
function of y, where y belongs to a suitable Banach space ft. We take 

h = {y = (y(q)) | y(q) e C 2d , ||y|| = E |y(?)| < °° }■ 

<? 

Let -B(r ) be the open ball of radius r in ft, centered at the origin and let H°°(B(r ), h) 
denote the Banach space of analytic functions || w : B(tq) — > ft, equipped with the 
supremum norm, which we shall denote by |||w|||. The bound implies that Wo G 
H°°(B(ro), ft) for ro small enough, but before stating this, it is convenient to encode the 

decay property of the kernels inherited from the estimate ( |3.1|) as a property of the 
functional wq- 

For that let Tg denote the translation by f3 G R d , (r J gF)(0) = F(0 — /3). On ft, is 
realized by (rpy)(q) = e l ^' q y{q). It induces a map w i— >■ wp from H°°(B(r ), ft) to itself if 
we set 

wp{y) = Tp(w(T_py)). 

On the kernels u/™**, this is given by 

wf\q ]qi ,...,q m ) = e i ^-^w^\q;q 1 ,...,q m ). (3.4) 
and makes sense also for (3 G C d . We have 

oo 

IIMH < E sup E e_ImM9 " E%) k m) (9;?i ) --- ) 9m;C)iC 

m=0 8l>— >9»» g 

Combining this with the bound (|3.1|) we can summarize the above discussion by 

Proposition 1. There exists tq > 0, a > and D < oo, such that wop G H°°(B(ro), ft) 
and it extends to an analytic function of (3 in the region |Im/3| < a with values in 
H°°(B(r Q ), ft) satisfying the bound 

\\\w p\\\ < D\X\. (3.5) 
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Moreover, wop is analytic in ( for \(\ < tq. 

Remarks. 1. The analyticity of wop in a strip in C d centered on R d followed from the 

exponential decay of the kernels w^^q; qi, . . . , q m ) as functions of q — J2qj- In order to 
show that the solution of Eq. ( |2.3| ) is real analytic in 0, we shall need to inductively estab- 
lish such a decay for the kernels of the Fourier transforms w n of effective functionals 
W n . This will follow once we establish the analyticity of w n p in f3 with uniform bounds. 

2. From now on, analyticity in ( will always be understood to hold for \(\ < r . 

We finish this section by collecting, for convenience, some standard properties of 
bounded analytic functions defined on open balls in Banach spaces (that are identical 
to those of analytic functions on finite dimensional spaces, see ||). Let h, h' , h" be 
Banach spaces, B(r) C h, B(r') C h! and w { G H°°(B(r), ti), w G H ca (B(r'), h"). Then 

Composition property. If |||u>j||| < r' then w o Wi G H°°(B(r), h") and 

\\\w o Wi\\\ < \\\w\\\. (3.6) 

Inequalities. First of all, one deduces from the Cauchy estimate that for ri < r', 

sup WDwixn^^Kir'-nr'WlwlW, (3.7) 

||a;||<ri 

where C(h'; h") denotes the space of bounded linear operators from h! to h". Taking 
r i — \ r ' ' i w e infer that if |||wj||| < \r' then 

\\\w o wi — w o w 2 \\\ < — \ \\w\\\ || \Wi — W2\\\ ■ (3.8) 
fe-i 

Moreover, if 5kw(x) = w(x) — J2 77 D e w(0)(x), then 

e=o 

sup ||<Wa;)|| < ^-|| k||| (3.9) 

||z||<7r' 

for < 7 < 1. 



4 Inductive bounds 



In this section, we first define the operators T n used in the RG transformations. Then, 
we state an easy result, namely that the fixed point equations ( 2.13 ) and ( p. 15 ) may be 
solved for any n, if we choose A small enough in an n-dependent way. Then, we define the 
spaces in which our RG equations fl2.13|) are eventually solved inductively for |A| small 
uniformly in n and state precisely our inductive assumptions. 



In the Fourier variables, the fixed point equations ( |2.13| ) and ( |2.15| ) may be written in 
the form 



w n p{y) = w {n _i )p (y + T n _ l w nP (y)). 
w n p{y) = w 0/3 (y + T <n w n(3 (y)) : 



(4.1) 
(4.2) 
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where we have introduced (5 assuming that the operators T n are diagonal in Fourier space 
and hence commute with rg. Similarly, the equations Q2.14j ) and (|2.16|) translate in the 
Fourier space to the relations 

fnp(y) = f(n-i)(}(y + r n _ity n/9 (j/)), (4.3) 
fnp{y) = y + T <n w nl3 (y). (4.4) 

To define the operators T n , we shall use an analytic partition of unity (x n ) dividing 
the positive line into scales. Let 

Xo (k) = 1 - e"* 6 , Xn(«) = e- (r?_n+lK)6 - e- (rr " K)6 for n > 1. (4.5) 

Clearly, 

oo 

EXn(«) = l, Xn(/0=Xi(7r n+ V> forn>l. (4.6) 

n=0 

Note that K _6 x n (/t) are entire functions of k and that, for |Imft| < B and £ = 0, . . . , 6, 

|«-^Xo(«)| < C, (4.7) 

< Ce"* |K|6 , (4.8) 

for some -B-dependent constant C. Define 

Tn(g, q') = Xn(u ■ q) G Q (q, q') = 7 n (w ■ q) 5 m >, (4.9) 
where the matrix 7„ is of the block form: 

7»(«) = Xn(«)K- 2 (_^ 'o) (4.10) 

and we denote the kernel of an operator a in C(h; h) by a(q,q') G End(C M ). We shall 
also need below more general operators T n (n) with shifted kernels, 

Fn(K)(q,q') = 7 n (w ■ q + k) 8 q>q , . (4.11) 

It follows easily from the bounds ([4.7]) and ( (4.8|) that for n > 1, 

lirn-i(«)|| £(hih) , l|r<„(«)|| £(h!h) < c r 2n if |«|<^s (4.12) 

with a new constant C (say, twice bigger). 

Our goal is to show that w n and /„ exist as analytic functionals provided that |A| is 
taken small in an n-independent way. For later purposes it will be useful to prove this first 
for |A| small in an n-dependent way. Although this is very easy to carry out, it illustrates 
an important part of the main analytic argument in the general step. 

Proposition 2. For any sufficiently small r > 0, |A| < A n (r) and |Im/3| < a, the 



equations Q4.2Q have a unique solution w n p G H°°(B(r n ), h) with 

\w nP \\\ < D\X\, (4.13) 
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where D is as in Proposition 1. The maps f n p defined by Eqs. ( |4.4[ ) belong to H°°(B(r n ), h). 
They satisfy the bounds \ \\f n p\\\ < 2r n . Moreover, w n and f n @ are analytic in X, [3 and 
C, and they satisfy the recursive relations (|4.1| ) and Q4.3|) , respectively. 

Postponing the proof to the end of the section, we shall state the bounds for w n that 
will be inductively established for |A| small in an n-independent way. Due to the smallness 
of [a; • g| in the n th scale, 7„ will have very different effects in the variables 9 and j in 
y = (9,j)- It will be therefore convenient to choose n-dependent norms for n > 1. Let us 
first do it for C 2d by defining 

|(*i,*2)L = \zi\ + ^\z 2 \. (4.14) 

We shall use the notation | ■ | n m for the matrix norms induced by viewing a 2d x 2d 
matrices as maps from C 2d with the norm | • |„ to C 2d with the norm | ■ | m . Next we set 

lblL = El^)Le^ nM . (4.15) 

The weight e v "l^' 9 ' will facilitate dealing with non-dangerous large denominators \u • q\. 
For w, it turns out to be useful to introduce the norms 

IML = £K9)Le-^l<^. (4.16) 

Let h± n denote the corresponding Banach spaces. Note the natural embeddings for n > 2 

— ► h, h — ► h- n+ \ — ► h_ n (4-17) 

with the norms bounded by 1: 

II ' II < II ■ L-l < II -L, II • ll-n < II • ll-n+1 < II • II • (4-18) 

For n > 2 (but not for n = 1), the operator T n _i or, more generally, operators r n _i(/«) 
may be considered as mapping h_ n into h n . Indeed, it follows easily with the use of bound 
that 



\\Tn-M\l n ., n < C V - 2n if l/^Tf-^B (4.19) 

with a new (n-independent) constant C. 

To simplify notations, we shall denote by B n the open ball in h n of radius r n and by 
A n the space H°°(B n , h_ n ) of analytic functions on B n with the supremum norm denoted 
by HI ■ |||. Finally, for a linear operator M : h n — > h m we use the abbreviated notation 
for the norm in £(h n , h m ). Due to the embeddings (|4.17|) , we may regard the maps 



Wnp, whose existence for sufficiently small |A| is claimed in Proposition 2, as belonging 
to A n . Note that both sides of relation (|4.1|) are well defined for such maps due to the 



bound ( |4.19|) and that their equality is implied by the results of Proposition 2. The next 
proposition states that, viewed as A n - valued functions of A, u>„^'s may be analytically 
extended to an n-independent disc |A| < Ao (provided we restrict somewhat the strip of 
(5). It also lists the properties of such extensions. 
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Proposition 3. (a) There exist positive constants r and \q with r < f] A such that, for 
|A| < A and |Im/?| < a n , where 

ai = a, a n = (1 — n _2 )a n _i, n>2, (4.20) 

there exist solutions w n p = w of Eqs. Q4.1|) belonging to A n , analytic in X, (3 and ( and 
coinciding with the solutions w n p of Proposition 2 for |A| < A n (r). 

(b) Writing 

w(y) =w{0) + Dw{0)y + 5 2 w{y), (4.21) 

we have 

\\Pw(0)\\-n < er 2n , (4.22) 

lilMil < er* B , (4.23) 



Dw(y)\\ n] _ n < er, 2n . (4.24) 



where e — > as A — > 0. 

(c) 



Remarks. 1. If we rescale the maps w n by introducing w n (y) = 7]~ 2n r~ n w n {r n y) then 
it follows from the above statements that w n p = are analytic maps from a unit ball in 
h n to h_ n and w(y) = ?D(0) + Dw(0)y + S2w(y) with 

||JM0)||_ < e V - 2n r n , \\\S 2 w\\\ < e V - 2n r? n } \\Dw(0)\\ n> _ n < e. 
Hence with growing n, Pw becomes an approximately linear map. 

2. Let us explain the idea of the proof of Proposition 3. Consider the linearization of 
Eq. §J: 

w n = m„_i + Dw n -iT n ^iw n + . . . (4.25) 
In order to solve the above equation one has to invert the operator 1 — Dty ri ,_ 1 r n _ 1 : 

W n = (1 - DWn-iTn^^Wn^i + ... 

However, operator r„_! is of order rj~ 2n as a map from /i_ n to h n (recall the bounds 
( f4.19| )) and we need to show that Dw; n „i is effectively of order rj 2n as a map from h n 
to h- n , which is, essentially, what Eq. ( f4.24| ) says with n shifted to n — 1. Altogether, 
Dwn-iTn-i remains of order e as a map from h_ n to /i_ n (this motivates also our choice 



of the norms) and ||(1 — Dw n -iT n ^x)\\ n ._ n < l + 0(e). In the proof of the estimate ( |4.24j ), 
we shall need the Ward identities discussed in Sect. 5. This is the only subtle part of our 
argument. Indeed, once the bound ( [4.24[ ) is shown, the rest of the proof of Proposition 



3 reduces to the standard Banach Fixed Point Theorem combined with the Diophantine 
property of u>. 
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The latter is used in the following way (which is similar to the way it enters the 
standard KAM proof): upon iteration, we consider smaller and smaller \u ■ g|'s, of order 
if 1 . This means \q\ is of order -r] n / T , by the Diophantine condition (|1.6|) . On the other 
hand, the introduction of the parameter (3 in ( |3.4j ) allows to preserve the exponential 

decay of the kernels Wo"'\q\ q\, . . . , q m ) in the size of \q — J2Qj\- By shrinking at each step 
the analyticity region in (3 we show that the leading contribution to w^s given by w n (0) 
(see Eq. ( |4.21| )) contracts for q ^ 0. Actually, ||Piu n (0)||_ n decays super-exponentially in 
n, see the estimate ( |6. 14| ) below, which explains why we can choose r as small as we want. 

Finally, the bound ( J4.23[ ) is easy to understand. By definition, 52W n and its first 
derivative vanish for y — 0, and the norm 1 1 l^^nl 1 1 is defined by taking the supremum 
over balls of radius r n , hence one expects IH^^nlH to be of order (1 + 0(e)) n r 2n by the 
Cauchy estimate ( |3.9| ) (the weaker bound fl4.23p is sufficient and is a convenient way to 
control the non-linear corrections to the iteration). Recall that, eventually, we construct 
our solution as a limit of X n = F n (0), for which we need to control w n (y) only for y = 0, 
see Eq. ( |2.16j ). Thus we can let the radius r n of the ball where our estimates hold tend 
to zero. 

3. Combining all the bounds, we get 

||K-(1-PK(0)||| < Ceif n . (4.26) 



The zero mode part (1 — P)w n (0) of w n (0) will be controlled later, see Eqs. (|5.3|) and the 
second of Eqs. (|7.4j ) below from which it follows that it is of the form (0, £ n ) where £ ra = 
O(X) converges in R d when n — > oo. Note that, since w n is multiplied by T n _! = T n _iP 
in the argument of w n _i in Eq. ( |4.1| ), the constant mode (1 — P)w n (0) may be decoupled 
from the iteration and we do not need to control it in order to prove Proposition 3. 

4. We choose the constants as follows, rj < 1 has been fixed first. B which enters the 
estimates ( [H2j ) and ( |4.19| ) is chosen then large enough (see Eq. ( |6.35p below). Given 



those, r and then Ao are chosen small enough. It should be emphasized that all quantities 
that are bounded by an n-dependent power of r are easy to estimate and that these 
estimates do not involve the Ward identities, the latter entering only in bounds with rj 211 . 
Finally, we denote by C a generic constant, independent of e and n, which may vary from 
place to place. 

Let us end this section with the easy 



Proof of Proposition 2. Consider the fixed point equation (4^) and write it as 
w = J-'iw) for w = w n/ 3 and 

J r (w)(y) = w op (y + V <n w{y)). 

Let B n denote the closed ball composed of w e H°°(B(r n ), h) with |||w||| < DX n (where 
D is as in Proposition 1). Choose A n so that Cr]~ 2n DX n < r n with C as in the bounds 
( f4.12|) . It follows from the latter that for w G B n and y G B{r n ) C h, 



\\y + T <n w(y) \\<r n + C V ~ 2n DX n < 2r n < \r (4.27) 

for r sufficiently small. Thus J^{w) is defined in B(r n ) and, by Proposition 1, ||jFu;(?/)|| < 
D\X\ < DX n . Hence T : B n — » B n . For tOj G B n , i = 1,2, use the property flOl) to 



12 



conclude that 



| H-Ftwi) - J r (w 2 )\ 1 1 < ^-Cr] 2n L»A n |||w;i - u> 2 ||| < 1 1 l^i - w 2 \\\ < ^\\\w 1 -w 2 \ 

' '0 * 



for r as in the estimate ( ^4.27| ), i.e. that T is a contraction. It follows that Eq. ( |2.15| ) 



possesses a unique solution w n p in B n satisfying the bound ( [4. 13| ) which, besides, is analytic 
in A, /3 and (. Consider now for n > 2 the map T'\ 

T\w)(y) = w p(y + T n ^iw n/3 (y) + r <n _iw(y)). (4.28) 

Again, T' is a contraction in B n since, for ||y|| < r n , we have ||y + r„_iti> n /3(y) + 
r <n _iw(y)|| < 3r n < |r for r sufficiently small. But Eqs. ( |4.2| ) imply that w n/ j and 
W(n-i)/3°(l+r n _iit; n/ g), both in B„, are its fixed points and, consequently, they have to coin- 
cide. Hence the recursions fl4.1|) follow. By virtue of the estimate ( |4.27| ), ||y + r <n u;(y)|| < 
2r n for y e B(r n ). By definition ( |4.4j ), this gives the claimed bound on |||/n/3|||- The re- 
cursion ( |4.3| ) follows easily from Eq. ( [4.1| ). □ 

5 Ward identities and cancellation of resonances 

The goal of this section is to prove the properties of the maps W n which will be essential 
in the proof of part (c) of Proposition 3, see Lemma 2 below, and in the proof Theorem 
1, see Lemma 1. These properties, which are proven by a simple integration by parts, 
result from the symmetries of Wq and will be encoded in the identities which, in the QFT 
formulation of the problem explained in Sect. 8, can be interpreted as the Ward identities 
corresponding to the translation symmetry. Their function in the proof is to guarantee 
a partial cancellation of the repeated resonances that plague the Lindstedt series, see 
Sect. 9. 



Indeed, as we shall see, the subtle part of the proof of the estimate (f4.24|) reduces to a 



bound on the diagonal elements of the kernel Dw n (0)(q, q) of the derivative Dw n evaluated 
at zero. Our strategy will be to show that this kernel is actually a function of u • q only 

and is of the form ( °^ w q ^ ^ °onf ) ior ou-q small. This, combined with our choice of the 

norms (|4.14j) and (}4.16|) , will then be used to imply the estimate ( f4.24|) . To show such a 



behavior of fw„(0)(g, g), we shall compute certain of its derivatives at lu ■ q = and show 
that they vanish. This is the role of the Ward identities proven here. Since q is in Z d , to 
make sense of such derivatives we need to introduce a smooth interpolation of Dw n (q, q) 
viewed as a function of u ■ q. This is the role of the functions 7T n defined after Lemma 1 
below. 

For simplicity, we shall first state and prove the Ward identities for the maps W n 
constructed in Proposition 2 for |A| < A n (by analyticity in A, they will also hold for W^s 
which will be constructed in Proposition 3 for |A| < Ao). The basic identity reads (recall 
that Y l = 6* for i < d) 

I W^<f>;Y)d4>= [ Y l (4>) d^W^frY) d<p iox i<d (5.1) 

JT d JT d 
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or, in the Fourier language, 

<{o-,y) = iY,y l te)j w ti-uy)> ( 5 - 2 ) 

where on the right hand sides the summations over the repeated index I from 1 to 2d are 
understood. Let us check first the n = case, see Eq. (|2.1|), 



J (<9,f/)((0, C) + Y((f>)) d<f> = Jdi [f/((0, C) + Y((f>))} d(f> 

\d l U){{(j )) C)+Y{ ( l ) ))d i Y l {4>) d<f>. 



The first term on the right hand side vanishes and the second one yields the claim by 
integration by parts. For n > 1, using the relation ( [2.15[ ), we obtain 



W^cf>; Y)d<P= I W*(cf>; Y + T <n W n (Y)) d<P 

(Y l + {T <n W n )\Y)){<P) d^Wkfo Y + T <n W n (Y)) d(j> 
Y\<f>) d^W^fr Y)d<P+ [ {T <n W n {Y))\ct>) d+iW^fr Y) 



The last integral can be written as / T<' n (0 — 0') W% (</>', Y) d^W^tp; Y) dtp d<p' and two 
integrations by parts and the symmetry of T <n show that it is equal to its opposite, hence 
that it vanishes. 

To derive the consequences of the Ward identity ( |5.2|) used later, evaluate it first at 
y = 0: 

Lemma 1. ^(0; y)\ y=0 = for all n and i < d. (5.3) 



The next identities involve Dw n . Let us first introduce smooth interpolations of the 
diagonal parts of the kernels of the derivatives Dw n p of the maps w n p constructed in 
Proposition 2 for |A| < A ra . These derivatives are given by the formula 

Dw n/3 (y) = [1 - Dw p(y n ) Y < n ] -1 Dw p(y n ) (5.4) 

with y n = y + T <n w n p{y). This is obtained from the ^/-derivative of Eq. ( [4.2|) . We shall 
now proceed to show that the kernel Dw n p(q,q'] y) on the diagonal q = q' depends on 
q only through uj ■ q. For this purpose, let us introduce the continuous automorphism 
t p : C(h; h) — > C(h; h), for p G Z d , shifting both arguments of the kernel of an operator a 
by p: 

(t p a)(q,q') = a(q + p,q' + p). (5.5) 

For n = 0, t p Dwo = Dwop, as follows from the explicit form ( |3.3| ) of the Taylor coefficients 
of Wq. Applying t p to Eq. (|5.4|) , we obtain the relation 

t p Dw n p(y) = [1 - Dwopiy^tpT^}' 1 Dwo/3(y n )- (5.6) 
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Note that t p T <n = T <n (uj -p), see the definition ( |4.11| ). It follows that t p Dw n/ 3 depends 



on p only through the scalar product uo-p. Denote Dw 0/ 3(y) = vr 0( g(?/) and define for n > 1 
and \k\ < i] n B, 

K n p{K\y) = [1 - ir 0f3 (y n ) T ^(k)}' 1 Tc op (y n ) . (5.7) 



Since, by the inequalities (|4.27p , \\y n \\ < |?"o for y e B{r n ) C h, Proposition 1 and the 
Cauchy estimate (ft.7|) imply that ||7To,3(l/n)|| £(h . h) < ^^|A|. It follows then easily that 
^n^i^'i y) is analytic for |A| < A n (r), |Im/?| < a, |«| < ^ n i? and y G B{r n ) C fa with the 
norm arbitrarily small if A — > 0, e.g. 

hnp{^y)\\ c{h ., h) < |A| 1/2 . (5.8) 

Comparing Eqs. ( |5.6|) and (|5.7|) , we infer that 

t p Dw n p(y) = 7i n/3 (uj ■ p; y) and t p ir n p(K) = 7i n/3 (K + u-p), (5.9) 

whenever defined, i.e. that / K n p{K) is a smooth interpolation of t p Dw n p. Note an ex- 
plicit expression, which we shall need later, for the ^-derivative of vr n( g(/t) obtained by 
differentiating Eq. (|5.7|): 

d K n nl3 (K;y) = ir nl3 (K;y)d K T <n (K)ii nl3 ( y K;y). (5.10) 

As is easy to check, the maps ir n satisfy for \\y\\ < r n and \k\ < r] n B the recursion relation 

7i nf3 (K;y) = [1 - n( n -i)p(K; y) r n _i(/t)] _1 7r (n _i )/3 (fi;; y), (5.11) 

where we have denoted y = y + T n _iw n p(y). 

The second consequence of the Ward identity is 
Lemma 2. ^'(0, 0; k; y)\ «=o =0 for i or j<d (5.12) 
3^(0,0; k; y)\ K =o =0 for i and j < d and t < 2. (5.13) 



y=0 



Remark. Eqs. (gl2|) , (|5l3D and imply that vr n (0, 0; k\ y)\ y=0 = ( o\1) )> the 



fact which will be used in an essential way in the next section. 

Proof. First, evaluating the derivative of Eq. ( |5.2| ) w.r.t. y at y = 0, we obtain: 

£>< z (0, g;0) = ig*u4(-g;0) for z < d. (5.14) 
The kernels of Dw n possess the symmetry property 

D<(g, <?'; y) = Dw^-g 1 , -q; y) (5.15) 

which follows by Eq. (|5.4j ) from the similar property of Dw l l (q, q'; y). The latter comes 
from the fact that Dwq is the symmetric second derivative of the functional A / U((<f), () + 
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Y(4>))d(p and Eq. ( |5.15| ) encodes, at least formally, the analogous property of W n , see 
Sect. 8. More generally, 



(q, q'\ «; y) = K. v) 



(5.16) 



where we denote the kernels of the maps ir n (n;y) by ir n (q, q'; k; y). Setting q = in 
Eq. (|5.14j) and using the symmetry (|5.15|) , we obtain the first claim 



D<(0,0; y)\ y=0 = < J (0,0; k; y)\«=o = for i or j < d 



y=0 



(5.17) 



For the second claim, we use the relation ( |5.10| ) for d K Tt n which, written in terms of the 
kernels n n (q, q'; k\ y), yields: 

d K K n (q, q'; k; y)\**o = V Dw n (q, q"; 0) <9 K 7 <n (u; • q") Dw n (q" , q'; 0) . 

v=° 



In particular, for q = q' = 
9 k < 3 '(0,0;k; y)\ K =o 



y=0 



E^<'(0, ?; 0) • ?) ^^'(0, -g; o) ; 



where we have also used the symmetry ( 5.1 5|) . Finally, for i,j < d, substituting the Ward 
identity (|5.14j) , we obtain the relation 

d K 7r»(0,0;K;y)\ K =o = w l n (-q\ 0) d Kl<n (u ■ q) ll 'w^(q; 0). (5.18) 



The right hand side of the last expression is symmetric in indices i, j but, by the symmetry 
( |5.16|) , the left hand side is antisymmetric, hence zero. We can also see this more directly 

since, as follows from Eq. (|4.10|) , d K ^ <n {n) = (^(2) o ) w hh a odd and b even. Thus 

the expression summed on the right hand side of Eq. ( |5.18 ) is odd in q and the sum 
vanishes. □ 



6 Proof of Proposition 3 



The proof contains two different parts. The inductive proofs of parts (a) and (b) are 
straightforward applications of the Banach Fixed Point Theorem and the contraction in 
n follows easily by combining analyticity with the Diophantine condition. The proof of 
part (c) may also be divided into two parts. First, one controls Dw n (y) — Dw n (0) and 
the off-diagonal elements of Dw n (0)(q, q'). This is also straightforward and similar to the 
proofs in parts (a) and (b). For Dw n (0)(q,q), we use Lemma 2 above. 



For n < n with fixed no, the bounds ( |4.22j ) and ( 4.23 ) as well as ( 4.24 ) follow imme- 
diately from the estimate ( 4.13]) and the inequalities | • || < || • || n , || • ||_ n < || ■ || by taking 
A small enough. For n > tlq with no large enough, we shall proceed inductively. It will 
be convenient to modify slightly the simplified notations of the text of Proposition 3 and 
so, below, w will stand for tf(„_i)/3 and w' for w n /s. Finally, T will stand for T n _i. 
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Proof of (a). Consider the recursive equation Q4.1| ) for w' and use the decomposition 
( [4.21|) to rewrite it as 

w'(y) = w{0) + Dw{0){y + Vw'(y)) + 5 2 w{y + Tw'{y)) 

from which we deduce that 

w'(y) = Hw(0) + HDw(0)y + u(y), (6.1) 

where 

u{y) = HS 2 w(y + Tw\y)) = HS 2 w(THw(0) + Hy + Tu{y)) (6.2) 
with H = {1- Dw(O)r)- 1 and H = 1 + THDw(0) = (1 - rDiu(O))- 1 . 

In the inductive step, first we assume that w sa tis fies th e bou nds ( |4.22j ), ( |4.23| ) and 
( 14.241) with n replaced by n - 1. The bounds (|4.24| ), ([4.19|) and ( [4.18|) imply then that 
the operators H and H are well defined with 

\\H\\ , \\H\\ , , < 1 + Ce < 2, (6.3) 

II II— n+1; — n+1 * II lln— l;n— 1 — 1 — * \ / 

for e (i.e. |A|) small enough. 



We solve Eq. (|6.2| ) using the Banach Fixed Point Theorem. Given its solution u, the 
existence of w' satisfying Eq. (|6.1|) follows. To solve Eq. (|6.2|) , we consider the map Q 
defined by 

g(u) (y) = H5 2 w{y) with y = THw{0) + Hy + Tu{y) . (6.4) 

We claim that Q is a contraction in the ball 

B = {u E H°°(Bs, fr-n+i) I HHII <2er^ n " 1) }, (6.5) 

where Bs C h n _i is the open ball of radius r n ~ s for < 5 < 1 and r < r(5). Indeed, 
for ||y|| n _ 1 < r n ~ s , the inequalities ([4.19 ), (|6.3|) and ( |4.22|) combined with the identity 
THw(0) = THPw(0) (HP = H follows from TP = T and the definition of H) and the 
inequalities (|4.18| ) imply that for u G B, 

||y|| n -i < 2Cr ] - 2n tr 2{n - l) + 2r n - s + 2Cr]- 2n er^ (n_1) < ir"- 1 (6.6) 

if r is small enough. Thus y G -B n -i, the domain of definition of w and hence of S 2 w. It 
follows that Q(u) for w in the ball (|6.5|) may be considered as an analytic map of y G B$ 
with values in h- n +i. Moreover 

||^(^)(y)||_ n+1 < 2|||MII < 2er»<"- 1 \ (6.7) 
where we used the bounds ( |4.23| ), ( [4.18| ) and (|6.3|) . Hence Q : £> — >• i3. 



17 



Tui(y) and Ui G B, i = 1,2. By inequality (|6.6| ), |||/i||_„ +1 < \r n 1 and so the bounds 



To prove that Q is a contraction, we use the estimate 
and Ui G B, i = 1, 2. By inequa 
((03|), (CT) and imply that 

\\\GM -GMWl = sup \\H8 2 w(m) - HS 2 w(y 2 )\\_ n+1 

< 4r- n+1 \\\5 2 w 

-' 2) ' ]U Ul -u 2 



for y i {y)=THw{0)+Hy- 



< 4er 2 



|(n-l) 



= sup \\H5 2 w(yi) 
yeB s 

sup ||^i -y 2 |L_i < 4er2 (n_1) sup 



2/2 1 



< -|||mi - u 2 | 



for r and e small proving the contractive property of Q on B. Hence the existence of the 
fixed point u G B of Q solving the equation ( |6.2| ) and thus of w' : B$ — > ^- n +i given by 
Eq. (|6.1|) follows. Using the natural embeddings ( [4.17|) , we may consider B n as a subset 
of Bs, and w' may be regarded as an element of the space A n of analytic functions from 
B n C h n to h_ n . Note also that, since y = y + Tw'(y) (see the argument of 5 2 w in 
Eq. (|6.2|)), the inequality Q6.6Q may be rewritten as 



\y + Tw' 



, < v- 1 

n-1 — 2 



for y E B s 



(6.8) 



which implies that y + r«/(?/) G B n -\ for such |/. 



It is easy to see inductively that for |A| small enough, the maps w n p constructed in 
Proposition 2 give rise via the decomposition (|6.1|) to it's which solve the fixed point 
equation (|6.2| ) and belong to B. It follows that w' constructed above coincide for small 
|A|, and hence for all A in the common domain of definition, with w n p of Proposition 2. 



(6.9) 
(6.10) 



Proof of (b). Using the decomposition ( |6 . 1|) , we may write 

w\y) = w'{0) + Dw'(0)y + S 2 w'(y), 

with 

w'(0) = Hw(0) + u(0) , Dw\0) = HDw(0) + Du{0) 
and 5 2 w' = 5 2 u. Let us first iterate the bound ( |4.22|) . Note that 

Pw'{0) = PHPw(0) + Pu(0) 
since H = HP. As u G B, see the definition (|6.5|) , we have, 

(n-1) 



||P«(0)||_„ +1 < |K0)||_ n+1 < 2er? 



which, with the use of the estimates 
that 



(6.11) 

and Q02] ) to bound \\H Pw(0)\\_ n+1 , implies 



\Pw'(0)\\_ n < ||Pw'(0)||_ n+1 < 2er 2 ( n " 1 ) + 2er^ 



(n-1) 



(6.12) 



This seems weaker than what is needed to iterate the estimate (|4.22j ) but, as we 
shall see, one actually needs much less. The crucial pointQ is that the bound fl6.12j ) 

4 This is the first of the two places where we gain by working with w n @ for complex /3. 
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holds for |Im/3| < ct„_i, while we have to establish the estimate ( |4.22| ) for w ' only for 
|Im/3| < a n = (1 — n _2 )a„_i. For such /3, we infer, using the estimate (|6.12|) , that, for 

\w'(q',0)\- n+1 e^ n - 1 - an ^ e-^ 1 ^ < ||Pw'(0)||_ n+1 < e (6.13) 

if in the middle term we take w' corresponding to the value of (3 shifted to (3' = (3 — 
^ (a,, i-a„) It follows that 

\\Pw'(q;0)\\_ n < e J2e' n ' 2an - M e-^-^'"^. (6.14) 

oo 

Since a ra _i > n (1 — n~ 2 ) a > 0, the sum on the right hand side is clearly bounded by 

n=2 

Cn 2 . However, we may extract from the sum factors that are super-exponentially small 
in n. Indeed, for \u ■ q\ < i]^ n , we may extract from the first exponential under the sum 

a factor e~ ^ v " n ) due to the Diophantine condition ( |1.6| ). On the other hand, for 

\uj • q\ > r? 7 ™, we may extract a factor e -0 * 77 7 ^ from the second exponential. Hence the 
inductive bound fl4.22| ) follows for n > hq, and n$ large enough. 

Let us now iterate the relation ( |4.23j ) for <5 2 w' equal to 82U (see Eq. (|BTT|)). Recall that 
\\u(y)\\_ n+1 < 2er* (n ~ 1) for < r n ~ 5 , see the definition Q. The estimate Q 

with k = 2 and 7 = r s and the bounds (|4.18| ) imply then that, for \\y\\ n < r n , 

ll«W(v)||_ < i^er^. (6.15) 

Taking 5 > | and r < r(5), we infer that \\S2w'(y) \\_ n < era". This completes the 
inductive proof of (b) . 

Proof of (c). We shall use the maps 7r n/3 introduced in Sect. 5, and related to Dw n by 
Eqs. ( |5.9D . Recall that ir n p was constructed as a map from B(r n ) C h to C(h,h), for 
|A| < A„. With the use of the embeddings (|4.17|) , they may be viewed as maps from 
B n C h n to C(h n ; h_ n ). As we shall see, they can be extended to |A| < A . Let us split 
7r n( g(/t; 0) into the diagonal and the off-diagonal parts: 

n n p(n]0) = a n (K) +p n (K), (6.16) 

where a n (q, q'; k) = n n p(q, q; k; 0) 5 q>q i. Let us denote by D n the disc {k G C \k\ < r) n B}. 

Lemma 3. The maps tt = 7r n/3 : D„ x B„ — > C{h n ; h_ n ) extend analytically to |A| < A . 
Their extensions satisfy the relations 

t p Dw(y) = ii{u) • p; y) and t p 7i(K;y) = tt(k + u ■ p; y) , (6-17) 

whenever defined, and depend analytically on k, y, (3 and (. They obey the bounds: 

IIM«;i/)IL- < ^\ (6.18) 
lk(«)IU_ < W n , (6.19) 

||p(«)L_n < er*» (6.20) 
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where 8itt(k] y) = tt(k] y) — 7t(k; 0). 



Obviously, the bound ( [4.24| ) follows by combining Eqs. ( |6.17| ) for p = 0, ( p,18| ), ( |6.19| ) and 
( |6.20| ). In particular, we may use estimate ( |4.24| ) as an inductive hypothesis in the proof 
of Lemma 3. Hence, we are left with 



Proof of Lemma 3. We use the same notation as in the proof of Proposition 3. Let us 
show that we may use recursively the relations ( [5 -11 ) to construct for |A| < Ao the maps 
7r„/3 satisfying the identities ( |6.17| ). First, differentiating Eq. ( f4.1| ), we obtain the relation 



Dw'(y) = \l-Dw{y)Y}- l Dw{y)., 



(6.21) 



where y = y + Tw'(y). The right hand side is well defined for y £ B$ C h n _i, since, by 
inequality (|6.8|), y £ £> n -i for such y's. The bound ( |4.19| ) and the inductive hypothesis 
( |4.24| ) imply that \\Dw(y)T\\_ n+v _ n+1 < Ce. Denoting n = 7T( n _ 1 ) ( g we then define n' = 7T n/ g 
by relation ( |5.11|) , i.e. by 



K'{K\y) = [1 - 7r(/c; y)T(n)] 1 n(n;y). 



(6.22) 



The relations (|6.17|) for it' may be inferred by applying the automorphism ( |5.5|) to 
Eqs. ( |6.21D and (|6.22|) . The inductive hypotheses imply that, for k £ -D„-i and y £ B$, 
tt'(k] y) is defined in £(h n _i; h_ n+ i) and is an analytic function of its arguments. It fol- 
lows by induction that it coincides for | A| < X n with the map 7r n/ 3 constructed before, see 
Eq. (|5.7|) , also satisfying the recursion (|5.11| ). Note that 



7r'(«;0) = 

where yo = Tw'(0). The bounds 



7T 



!«;2A))r(«)] n(K;y 



|7T [K 



n — 1; —n + 1 



n— 1; — n+1 



< e ^(n-l) ; 



(6.23) 

(6.24) 
(6.25) 



follow easily: the stronger factor r 7< - n_1 ^ instead of the weaker one ^ 2 ( n_1 ) in the estimate 
of 5i7r' is supplied by the inductive bound (|6.18| ) for the difference 7t(k; y) — tt(k,; y ) = 
5iti(k,; y) — 5iir(n] y ). The inequality (|6.25|) is enough to iterate the bound ( |6.18|) . Indeed, 

the estimate (|3.9|) with k = 1 and 7 = r s permits to extract the additional factor 
and to obtain the improved bound (|6.18|) (if 8 > % and r < r(5)) for H^iTt^k; y)||„._„ for 
y restricted to B n C h n . 



Let us turn to the inductive proof of the bound ( |6.20|) for ||p / (K)|| n ._ n . We gain a very 
small factor from the restriction of the analyticity strip in /3, as before in the control of 
Pw'(0) in the proof of (b), see the inequalities (|6.13| ). First, by the definition Q4.16j ) of 
the norms, we have the following estimate of the kernel p'(q, q'; k): 

1^(9, «)L_x ; _ +1 e^-^^l^'l e -^ +1 fl-^+l-^l) < ||7r'(K;0)|L 1; _ n+1 < e, 



where in the middle term we have shifted (3 to (3' = (5 



(q — q'). The last 

inequality follows from the bound ( |6.24| ). Moreover, for an operator with kernel a(q,q'), 



(«n-l-"ti) 

k-5'l 
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|| a IL-„ < supX) \a(q 1 q')\ n ._ n e v "(l^l+I^VI) . This implies the estimate 

qi q 

The expression on the right hand side may be bounded by a factor super-exponentially 
small in n, which is much more than is needed. Indeed, we may extract from it an extra 

factor e~°^ ^ n ' for \lo ■ q\ and \u) • q'\ < v^ n and hence 7^ \u -{q — q') | < 2?7 7?1 using the 

Diophantine condition (|1.6|) and an extra factor e~ ^ * ) if \u • q\ > rf n or \u ■ q'\ > rf n . 
Hence the bound ( |6.20|) for the off-diagonal operator p'{n). 

We are left with the proof of the estimate ( |6.19|) for the diagonal operator (j'(k). Let 
us define 

s (z) = ^-^{0,0; Br]" 1 - 1 z)u n -\ (6.26) 
where u = (0 ° ) is a block matrix. Similarly, we introduce the matrix s'(z) related to 

a'. With the use of symmetry (|5l6l) , we write s\z) = ^ 2(n+1) (^2) )• We shall 

prove that, for \z\ < 1, 

\p' t (z)-p>z 2 -*\ < A\z\^ (6.27) 

with |^| < (1 - and A < ^, assuming inductively similar bounds for s(z). Note 
that such inductive assumptions, together with the identity |M| n _ 1; _ n+1 = |-u n_1 MM n_1 | 0;0 
for the matrix norms imply, in particular, the estimate 

|a(0,0, «)]„_!;_„+! < \erf n (6.28) 

for k G -D,t-i- The bound ( |6.27|) will follow from Lemma 2 expressing the cancellations of 
resonances. The leading Taylor coefficients pi of pi(z) are marginal in the RG terminology 
and the higher ones are irrelevant. The presence of lower order relevant Taylor coefficients 
would spoil the iterative bounds. They are, however, forbidden by the Ward identities. 
Let us pass to the details. 

Let us first prove the estimate ( |6.19| ) for a' assuming the bound ( |6.28| ). We shall split 

a'{K) = (t' q (k) +(t[(k) with o£(k) = [1 - <r(«)r(K)] -1 «7(K), (6.29) 

compare with Eq. (|6.23| ). Note that <j(k)T(k) is an operator diagonal in Fourier space 
and hence, so is ct' q (k). Since, by the inductive hypotheses ( |6.20|) and (|6.18| ), 

lk(K)-7r(K;y Q )|L_ 1; _„ +1 = \\p(K) + 8 1 Tr(K;y )\\ n _ 1; _ n+1 < 2er^, (6.30) 
it follows that 

K(*)IL-» < W\ (6-31) 
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for r small enough. We pass to the estimation of ct' q {k). Note that the bound QL19| ) 
together with the inequalities fl4.18|) and the inductive hypothesis ( |6.19|) imply that 



\T{k)<t(k) \\_ n ._ n < \Ct so that 

K(«)IL_» < 2|K«)|| ni _ n . (6.32) 
For operators a diagonal in Fourier transform, ||a|| n ._„ = sup \a(q, q)\ n ._ n e~ 2v n l w '^l. Hence 
it follows from the bound ( |6.32 ) that 



IKo(«)IU-n < SU P 2|a(g,g; K )L ; _ n e- 2 ^l^l. (6.33) 
<? 

For g with |cj • q\ < (1 — r))r) n ~ 1 B, we use for k G D n the equality cr(g, g; k) = <r(0, 0; k) 
with k = k + uj ■ q which follows from the second identity (|6.17|) (observe that for such g's 
and for k G D n , k G -D n -i)- By virtue of the inequality (|6.28|) , 

|<T(0,0;£)| n; _ n < \a(0, 0;«)|„_ 1; _ n+1 < ier/ 2 ". (6.34) 

Hence, for g with |u; • g| < (1 — f])r] n ^ 1 B, we may bound the expression on the right hand 
side in the estimate ( p.33| ) by jer] 2n . For \uj • q\ > (1 — r])7] n ^ l B 1 we instead extract an 
extra factor estimating 

2k(?,^)L ; -„e- 2 ^ N < 2||a( / ,)L_ 1 _ +1 e- 2 ^ 1 ^) 2B 

< e -2l-Hl-r,)*B erj 2(n-l) < l_ erj 2n (g gg) 

for 5 sufficiently large (this is the only place where B large is needed). Putting these 
estimates together with the inequality ( |6.31|) for a[, we infer the bound (|6.19| ) for ct'(k). 



We still have to iterate the crucial estimates (|6.27|) which is the only place in the proof 



of Proposition 3 wh ere we use the Ward identities. Writing Eqs. ( |6.29| ) in terms of s, see 
the definition (|6.26|) , we obtain 



s\z) = [l-(Cs)(zmzT\Cs)(z) + s[(z) 



with the "linearized RG map" C, 

(Cs)(z) = us(r]z)u, 

and 7(z) = vT n ^ n „ x {Bri n z)vr n = r)- 2n (Bz)- 2 Xx{Brfz)( _£ Bz lB z ), see Eqs. (£L|) and 



( |4.10|) . The estimate ( |6.30D implies that the remainder s[ satisfies the bound |s'i(z)| 00 < 



Cer2 n . Combining the definition Q4.5[ ) (which implies that Xi( z ) is °f order \z\ 6 for small 
z) and the inductive bound for s, we infer that \(Cs)(z) ! y(z)\ 00 < Ce\z\ A . Thus altogether 

\s'{z)-{Cs){z)\ 0i0 < Ce 2 V 2n \z\ A + Cer? n . 

The map £ preserves pi and contracts the constant A to r]A. The Ward identity, Lemma 
2, implies that = for j < 2 — i. Since \p[ — pi\ < Cei]~ 2n r^ n < and 

t] A + Ce 2 + Cer]~ 2n r^ n < A for r and e small, we infer that s' satisfies the bound ( |6.27| ). 
This finishes the proof of Lemma 3 and of Proposition 3. □ 
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7 Proof of Theorem 1 



We shall first show that X n = F n (0) = T <n W n (0) converges to a real analytic function X 
with zero average as n — > oo and that X solves Eq. (12. 



Recall that in Proposition 2, we have constructed for |A| < A n the analytic maps f n p 
from B(r n ) C h into h, satisfying the relations (|4.4j ) and ( f4.3|) and the bound |||/n^||| < 
2r n . They may be also viewed as analytic maps from B(r n ) C h n to /i. As such, they may 
be analytically extended to |A| < A for n > n by iterated use of Eq. (|4.3| ) if we recall the 
bound (|6.8|) . The extensions are clearly uniformly bounded (e.g. by 2r n °). Let us prove 
the convergence in h of x n p = / n/ g(0) obtained this way. The recursion (|4.3|) implies that 

Xn? = X{n-T)P + 5l/(n-l)/3(r n _iW nj 3(0)). (7.1) 

Using the estimate ( |3.9|) for k — 1, we deduce from the inequalities ( |4.19|) and (|4.22j) that 
the || • || norm of the 2 nd term on the right hand side of Eq. (|7.1| ) is bounded by Cer]~ 2n r 2n . 
Hence the convergence 

x n /3 — > xp in h (7.2) 

together with the bound ||xg|| < Ce uniform in the strip \lm/3\ < a' = «n™ =2 (l — n~ 2 ). 
The latter implies the pointwise estimate 

\x{q)\ < Cee~ a ' lql (7.3) 

and, consequently, the real analyticity of the Fourier transform X of x. 

For |A| < A n , Eqs. flU) and (|J) imply that 

x n = / n (0) = T <n w (x n ) and w (x n ) = w n (0). (7.4) 

In particular, it follows from the first of these equations that x n (q)\ q=0 = and from the 
second one and Lemma 1 (i.e. the Ward identities) that w l Q (q; x n )\ = for i < d. By 
analyticity in A, these relations have to hold for |A| < Ao- Since wq is analytic, we can 
take the n — > oo limit of Eqs. (|7.4|) , and infer that 

x(0) =0, x(q) =G w (q;x) for q ^ 0, and w (0; x) = (0, £). (7.5) 



The first 2 of these equations are the Fourier transformed version of Eq. ( |2.3j ). The 
solution x is an analytic function of A for |A| < Ao and it vanishes for A = 0. Recall that 
( is a parameter in wq, and thus x is also analytic in ( for \(\ < r . 



We still have to solve Eq. (|J) for C In view of the 3 rd of Eqs. Q and Eq. ([D]) it 
reduces to the equalities 

AtC + A / d I U((4>,C)+X( ( f)))d ( f) = or f d I U(((f),()+X(<f)))d<f) = (7.6) 

in the, respectively, non-isochronous and isochronous cases. We shall solve the above 
equations for ( by the Implicit Function Theorem. Note that A = and ( = satisfies 
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the non-isochronous equation and that the ^-derivative of its left hand side at these points 
is /i, which is assumed to be invertible. Similarly, A = and ( = solves the isochronous 
equation (since X vanishes for A = and we have assumed that / diU(4>, 0) dip = 0) and 
the ^-derivative of its left hand side at these points is / d]U(<p,0) dtp, which is assumed 
to be invertible. The existence of the local solution £(A) analytic for |A| < Ao with Ao 
small enough and vanishing at A = follows in the both cases. The resulting solution 
Z = X + (0, C) of Eq. ( |1.3| ) depends analytically on A for |A| < A and vanishes for 
A = 0. Its uniqueness up to translations ( |1.5|) follows from the fact that the equations 
( |2.3|) and fl2.4|) completely determine the coefficients of the Taylor expansion of their 
solution in powers of A (i.e. the Lindstedt series discussed in Sect. 9). This ends the proof 
of Theorem 1. □ 



8 Field theory interpretation 

Although arising in a problem of classical mechanics of d degrees of freedom, Eq. ( |1.3| ) 
has a field theory interpretation, see JTTJ, Let us consider the (action) functional S(Z) 



of maps Z = (6, J) : T d ->• R d x R rf 



S(Z) 



\j{4>) ■ fiJ(<P) - J(<P) ■ (u ■ fy)6(0) + XU(<P + 6(0), J(<P)) 



Note the translation symmetry 

SfaZ - (P,0)) = S(Z), (8.2) 

where, as before, (rgZ)(0) = Z(<p — (3) for (3 G R d . As is easily seen, our basic equation 
( |1.3|) coincides with the equation 5S(Z) = for the extrema of functional S. The map Z 
with / G(0) dip = minimizing S may, in turn, be interpreted as the limit H — > of the 
formal functional integral 

m >* = fei^ DZ, <8 - 3) 

where DZ denotes the formal Lebesgue measure on LQ(T d ,H d ) /Z d x L 2 (T d , R d ) where 
L\ is composed of maps with zero average. Indeed, in this limit, the integral should be 
localized at the minimum of the functional (B-ll). 



While Eq. ( |8.3|) has a purely formal meaning, we may gain intuition from it by some 
further manipulations. First, write Z = (0, () + Y with Y of zero average. Then 

= j[(o,o+^)] ^ Vo[Y '° ^ocn^(c) f84 , 

where 

V (Y, C) = A / U((<f>, C) + (8.5) 
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and we have used the quadratic part of the action functional to define the oscillatory 



•'measures" 



dji(C) =(det^e^ c d( (8.6) 
and dfiQ. The latter is formally given by 

d^{Y) = e-^^DY J j \-^ Y ' G ^ DY , (8.7) 

i.e. it is the "Gaussian measure" with mean zero and covariance iHGq on the space of 
maps Y with zero average. 

In quantum field theory, the idea of the renormalization group (RG) is to calculate 
functional integrals inductively. Let us explain a concrete realization of this idea. Let 

Gq = G\ + T 

be a decomposition of the operator Go into a sum of two operators. We may then write 
the measure dfiQ as a product measure according to the formula: 

J F(Y) dfi (Y) = Jf(Y + Y) dfnOO duo(Y), (8.8) 



where d\i\ is the measure given by Eq. (|8.7| ) with Go replaced by G\ and du is the 
Gaussian measure with mean zero and covariance ihT Q , both on the space of maps with 
zero average. Using this identity, we may rewrite Eq. (|8.4j) as 



^ = Je^^flWO <8 - 9) 



if we set 



VL(y,0 = [ e i Vo( - Y+Y <Q dv (Y) , (8.10) 



f (y c- 6) J[ y w+ywi^ y+?i ° d _^n (8U ) 

Fl{YXA) ~ — d M Y) ■ ( } 

This is the first step of the iterative procedure. After the n subsequent decompositions, 

n-1 

Co = G n + Tfc, 

one arrives at the expression 

7IM \ J [(0. + F n (Y,C, 4>)} ej V "' V '° dfeWjMQ 
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with 



v ^ = [ e* v ~-& + *'Q dv^Y) , (8.13) 



F " (y ' C;0) = /e ^™^ l( y) • (8 - 14) 

Such a procedure may lead to the calculation of the expectation (Z(<f>)) if we have a 
good control of the asymptotic behavior of the effective interactions V n and of the effective 
insertions F n . Note that the translation symmetry ( |8.2j ) goes through the RG iteration 



V n {r p Y - (/5,0)) = V n (Y) (8.15) 

if we consider in the definitions (|8.10|) and ( |S.13| ) fields Y with arbitrary averages and 
if the covariances T n commute with the translations Tp. The latter property has been 
guaranteed by the explicit construction of the operators T n (see Eq. Q4.9|) ). 



The inductive RG scheme described in Sect. 2 may be obtained as the formal h — ► 
limit of the above iterative calculation of the functional integral (|8.3| ). After the first step, 
we obtain 

V 1 (Y) = -\{Y Q ,Y Q l Y Q ) + V Q {Y + Y ), F X (Y) = Y + Y , (8.16) 
where Yq minimizes the right hand side of the first equation. Denoting 

we may rewrite Eqs. (|8.16|) in the differential form: 

Wx(Y) = W (Y + Y ), Y Q = T W (Y + Y ) (8.18) 

or, more conveniently, as the fixed point equation 

W,(Y) = W (Y + T W 1 (Y)) (8.19) 

whose solution determines F\\ 

Fl {Y) = Y + Y = Y + ToW^Y). (8.20) 

These are Eqs. ( |2.9| ) and (|2.10|) , respectively. Similarly, after n inductive steps, W n and 
F n are determined by relations ( |2.13| ) and ( p,14| ), the h —>■ versions of Eqs. ( |8.13| ) and 



(|8.14|) , respectively. The cumulative expressions ( |2 . 1 5| ) and (|2.16| ) may be obtained as 



those for W\ and F\ by replacing T by T <n = Tfc, i-e. by performing n RG steps at 

k<n 

once. 

The Ward identity ( |5.1|) of Sect. 5, which assured the partial cancellation of the re- 
peated resonances or, in the field theory language, the absence of the terms proportional 
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to / |6| 2 , / O • (uj ■ d)Q, J J • 6 (and to / 0) in th e effective interactions V n , is the 
infinitesimal version of the translation symmetry ( |3.15|) . It is obtained from the latter by 
the differentiation w.r.t. (3 % at (3 = 0, which yields 



svjx) 



E 



0. 



The integration by parts in the second term and the definition ( 8 . 1 7|) give then the identity 



9 Renormalization Group and Lindstedt series. 

In this section, we sketch the connection of our approach to the Lindstedt series and to the 
resummation of the latter by Eliasson. The Lindstedt series has a graphical representation 
and our RG method amounts to resumming at each step a particular subset of graphs. 
Let us first introduce the graphical representation. It will be convenient to use a common 
symbol Q for the momentum variable q and the vector index i, Q = (q, i) G Z d x 

{1, . . . , 2c?} = J and to write x(Q) = x % (q) and WQ m \Q . . . , Q m ) for the matrix elements 

of the kernels ti>o (<Zo; <Zi • • • , q m \ C) introduced in Eq. (|3.3|). Using the relations ( |2.3| ) and 
(|3.3|), we obtain 

oo m 

<Q) = E E* G o(Q,Qo) 4 m) (Q)II*(<W, (9-i) 

m=0 Q £=1 

where Go(Q,Qo) = Go(q, qo) u ° and Q = (Qo, . . .Q m ). The sum YT means that we sum 
over qi ^ (this is due to the projector P in Eq. (|2.3j )). Since is proportional to A, 
this formula yields a power series solution in powers of A of x(Q) obtained by regarding 
the equality ( |9.1|) as a fixed point equation solved iteratively, starting with x — 0. 

The resulting series can be conveniently expressed as a sum over tree graphs whose 
weights are as follows. Let T^ +1 denote the set of connected tree graphs T on m + 1 + k 
vertices v G {0, . . . ,m+k} = V(T) such that the first m + 1 ones, the "external vertices", 
have coordination number one. We shall assume that m+1, k > 1 and will call the first of 
the external vertices the root of T. The lines i of T are pairs i = (v, v'), v, v' G V(T) which 
we order assuming that the unique path going from the root to v' goes through v. The 
set of lines is denoted by £(T), lines t containing an external vertex are called external, 
£ G £e(T), and the remaining ones internal, £ G Ci(T). Given a function G : J 2 — > C 
and a collection (w( m ') m >o = w of functions : J' m+1 — » C symmetric in the last 

m-variables, we define the "amplitude" A(T, G, w, Q) of T. For this purpose, we assign 
variables Qi G J to the external vertices < % < m of T and variables G J to each 
internal line i and a vertex t> contained in it. We write = (Pe v ,Pev') f° r ^ = { v -> v ')- 
For m < v < m + fc, we set R v = Pg v if there exists an internal line (v 1 , v) or R v = Q 
otherwise (i.e. when (0, v) is an external line), and P„ = {Pi v \ £ G Cj(T), £ = (v, v')}, 
Qv = {Qi \ £—(v,i)e C E (T)}. We define: 

m+k 

A{T,G,w,Q) = ^tE* IT G ( p e) II ™ v l w^R*, P v , Q„), (9.2) 

P l&d(T) v=m+l 
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where Q = (Qo, ■ ■ ■ , Q m ) , the sum J2* runs over all pi v ^ 0, and m v + 1 is the coordination 

p 

number of vertex v. 

With this notation, the iterative solution of Eq. ( |9.1|) is given by 

x{Q) = E* E G o(Q, Qo) A(T, G , w , Qo) , (9.3) 
Qo TeTi 

where the sum is over all tree graphs with one external line, T\ = Ufc7f\ Comparing with 
Eq. ( |9.1|) , we also infer that 

w (Qo; x) = E A (T, Go, w , Qo) (9.4) 
TeTi 



Remark. A formal way to derive the identity (|9.3| ) is to start from the field theory 
formula QS.4|) , expand it in Feynman diagrams, and let h~ — > 0. In that limit, only tree 
diagrams remain (each line has a power of h~, while each vertex carries a factor h~ l ; so, 
all graphs except tree graphs are multiplied by some positive power of h), and we obtain 
Eq. 



As is well known since Poincare |20| , the series in Eq. ( |9.3| ) does not absolutely converge 
for nontrivial potentials U in the Hamiltonian (|1.1| ), i.e. StgTi \A(T, Go, wq, Qq)\ = oo 
(see e.g. for a simple proof). This is due to the presence of repeated resonances, namely 
of long sequences of lines (connected by vertices with coordination number two) many of 
which have the same small denominators to ■ pg. However, many trees contribute to the 
same order in A and, as shown by Eliasson 0, cancellations occur when one regroups 
terms in a suitable way. 

Let us now explain the renormalization group in the graph language. The key idea is 
a combinatorial identity which performs a partial resummation in Eq. ( |9.3|) : 

Proposition 4. Let G = G' + T. Then, 

E A(T,G,w,Q)= E A{T',G',w',Q) (9.5) 



with 



(P)= E A(T,T,w,P) (9.6) 

TeT m+1 



The proof of the identity (|9.5|) is quite simple. Insert G = G' + Y in the expression 
of A on the left hand side of ( |9.5| ) and decompose the tree T into a family of connected 
subtrees {T a } containing only T-lines and joined together by G"-lines. Let T" denote the 
tree which is obtained from T by contracting each tree T a to a point. Now, let us fix T' 
and sum over the {T a }'s; this leads to the identity (|9.5|) with the "renormalized" vertices 
w' given by Eq. (|9.6j ). 
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We may apply the identity (|9.6|), starting with w = wq, G = Go = G\ + r , w' = w\, 
and then inductively to w = w n -i, G = G n _i = G n + T n -i, w' = w n . It is easy to check 
that the w' n s so defined coincide with those constructed through Eq. (|2.13|) . Indeed, by 
the definition (|2.9|) , 

W X {Y) = W (Y + T W 1 (Y)). 
Writing both sides in Fourier transform and expanding in a Taylor series, we obtain: 

oo m 

MQ;y) = EE^ m) (^ p )Il^) 

m=0 P e=i 

oo m 

= EE^fftPJllfe^ + ^if!/)^)) (9-7) 



u 

m=0 P 



Now expand the product over £ in the right hand side replacing Wi(Q; y) in T wi(y)(Pi) = 
J2q ^o(Pi, Q) Wi(Q] y) by the right hand side of Eq. (|9.7|) . This leads, upon iteration, to 
a sum over trees with a number of external lines to which a factor y{P() is attached. The 
trees have an amplitude A(T, To, wq, Q), Q = {Q,P) and we can rewrite the right hand 
side of Eq. (|9.7|) , using the definition (plf ), as 



EE^ m (Q,p)II^)- 

m=0 P 1=1 

Comparison with the left hand side of Eq. ( |9.7|) shows that wi = w' defined here coincides 
with wi defined in Sect. 2. The same arguments apply inductively to w n . 



Note that the transformation (|9.6|) is rather easy to control (in the sense that, if w is 
small, w' is also small in a suitable norm) and could be used to give an alternative proof 
of Proposition 3. Indeed, and this is the main difference between our approach and the 
one of Eliasson, each transformation ( |9.6|) involves small denominators (in T n ) on only 
one scale. To see intuitively how to use this fact, consider first all the trees having only 
vertices with coordination number different from 2. It is easy to see that these trees have 
a number of vertices with coordination number one which is proportional to their total 
number of lines. Thus, one can control for those trees the small denominators on the 
lines, all of the same size, by the exponential decay in \q\ of w^(Q) and the Diophantine 
condition ( |1.6|) . The next observation is that we may reduce ourselves to those trees by 
resumming the contributions coming from the vertices with m = 1 (i.e. coordination 
number 2). This is where the problem of repeated resonances appears in this formalism. 
We obtain a series of the form: 



oo 



^(r nW «) fc r n = (i - r n w^y^n = HT n 

k=0 

which has to be controlled using the Ward identities, as in the proof of Proposition 3. 
Here we see again that this is the subtle point of the proof. 

Finally, let us translate in the tree language some of the formulas introduced in Sect. 2. 
Applying the resummation (|9.5|), (|9.6| ) to Eq. ( |9.3| ), we obtain: 

<Q) = E* E G o(Q, Qo) A(T, G u Wl , Q ) (9.8) 

Qo TeTi 
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and, inductively, 

x{Q) = E* E G o(Q, Qo) A(T, G n , w n , Q ) (9.9) 
Qo reTi 

which, upon the substitution Go = G n + T <n , may be rewritten as 

<Q) = E* E G n (Q 1 Q )A(T,G n ,w n} Q ) 
Qo TeTi 

+ E* E r <n(Q, Qo) A(T, G n , w n , Qo) . (9.10) 
Qo TeTi 

This corresponds to the decomposition 

X = F n (Y) = Y + T <n W n (Y), 

see Eq. ( p. 161 ), with Y satisfying Y = G n PW n (Y). Indeed, the solution of the latter 
equation may be written, in the same way as the solution ( |9.3|) of Eq. (|2.3|) , as a sum 
over the tree graphs. This gives the first term of on the right hand side of Eq. (|9.10| ). 
Moreover (compare with Eq. ( |9.4|) ), 

w n {Qo;y)= E MT,G n ,w n ,Q ). 

This gives the second term on the right hand side of Eq. ( |9.10|) . 

Note also that X n , defined by the equality X n = T <n Wo(X n ), see Eq. (|2 . 1 7|) , satisfies 
the equality 

Xn(Q) = E* E T <n(Qi Qo) A(T, T <n , Wo, Qo), (9.11) 
Qo TeTi 

which may be derived from the equation X n = T <n Wo(X n ) in the same way as we derived 
the identity (|9.3| ). Applying the relations (|9.5| ), (|9.6| ) inductively, we we infer that 

Xn(Q) =J2* T <n(Q,Qo)w { °\Qo) (9.12) 
Qo 

which is equal to the term in the second sum of Eq. fl9.10| ) corresponding to the tree with 
only one vertex. It is also easy to see why the other terms in Eq. ( |9.10| ) are small and 
therefore why shown in the proof of Theorem 1: in all the other terms, the 

w n ° (p)' s attached to vertices with coordination number one are multiplied by some G n , 
which forces the corresponding p to be large (using the support properties of G n , the 
Diophantine property (JO]) and the fact that the sums here run over p ^ 0). Then the 
exponential decay of w^Jp) implies that the contribution of those trees is small. 

Acknowledgments. We thank L.H. Eliasson, G. Gallavotti and V. Mastropietro for 
useful discussions. 
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